Monte Carlo calculation of the linear resistance of a three dimensional lattice 
Superconductor model in the London limit 
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We have studied the linear resistance of a three dimensional lattice Superconductor model in the 
London limit London lattice model by Monte Carlo simulation of the vortex loop dynamics. We find 
excellent finite size scaling at the phase transition. We determine the dynamical exponent z — 1.51 
for the isotropic London lattice model. 
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The fluctuation regime in high T c superconductors 
(HTCS) is expected to be sufficiently wide that critical 
fluctuations are observable fj],§]. In particular the con- 
ductivity is supposed to scale as a oc ^ 2 ~ d + z where 
£ is the correlation length and d is the dimension of the 
system. This scaling relation has been applied in recent 
experiments on YBCO in zero magnetic field ||. From 
which the value z « 2.6 and v ~ 1.2 (y is the corre- 
lation length exponent) was extracted. Accordingly an 
accurate determination of z and v in models of high T c 
superconductors is of great interest. The phenomenology 
of superconductors is described by the Ginzburg-Landau 
(GL) model. The model is to complicated to allow all de- 
grees of freedom to be included in calculations. Among 
the standard approximations of the GL model one can 
mention: the XY Villain and the lattice su- 

perconductor model in the London limit [p[-p6[ . 

In the present paper we determine z in the zero field 
London lattice model (LLM). The exponent z is known 
to be close to 3/2 in the 3 dimensional Ay-model, cor- 
responding to model (E) |17| ]. 

It is of interest to know whether the London model in 
which the spin wave degrees of freedom are integrated 
out is characterised by the same exponent. Equilibrium 
properties of the XY and the LLM for A = oo are known 
to be the same since they are connected through the Vil- 
lain duality transformation |J. However, the dynamical 
properties might not be the same. It is seen in other sys- 
tems where the spin degrees of freedom have an effect on 
the dynamics of the topological defects |fl8|| . However, as 
we show below, in fact the LLM has z — 1.5. This result 
is reassuring given that the model is used to study the 
dynamics of vortex systems in the relation to HTCS jl9| . 

Since the magnetic field H mag — we can limit our 
study to the isotropic system. We derive an expression 
for the resistance R, based on the Nyquist formula pc| ] 
for voltage fluctuations. From the Nyquist formula we 
derive a simple finite size scaling relation for the resis- 
tivity at the critical temperature T c and determine the 
critical dynamical exponent z. 

The LLM describes the vortex loop fluctuations of 
a bulk superconductor. The model originates from a 
Ginzburg - Landau description with no amplitude fluctu- 



ations and the spin waves integrated out within a Villain 
approximation. On a cubic lattice a vortex loop consists 
of four line elements forming a closed loop. 

The LLM is defined by the partition function Zona 
cubic lattice of side length L using periodic boundary 
conditions: 



Z = Trexph/lff] 
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where H is the Hamiltonian, the link variables q a rep- 
resent the vortex line elements. The are three kinds of 
q a , one for each direction e x ,e y and e z . The positions, 
of q a , are given by r-. The link variables q a i e { — 1, 0, 1} 
take are and ±1. The sum of q a over a unit qube 
equals zero. This is achieved by the trial updating al- 
gorithm, only to add closed vortex loops to the system. 
The Greens functions G a (r) |l3| are given by: 
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where k are the reciprocal lattice vectors, k x ,k y , and 
2irn/L, n = 0, . . . ,L — 1, k 2 — k 2 x + k 2 „ + n 2 



v 



k z 

and k x = sm{k x /2d), d is the side length of the unit 
cell and is set to rf = 1. The X x and A z are the 
bare magnetic penetration lengths in the x and z direc- 
tions. The coupling constants J x and J z determine the 
anisotropy of the model and are related to the screen- 
ing length by J z /J x = A 2 , /A 2 . In the work presented 
in this letter the penetration length is taken to be infi- 
nite, X x = X z = oo, we further restrict the model to the 
isotropic case J x — J z = 1. 

We simulate the model defined by eq. pi by the stan- 
dard Metropolis Monte Carlo method plj . The trial 
move consists of adding a closed vortex loop formed out 
of 4 link variables q. The loop is placed at a randomly 
chosen position and with one of the 6 different orienta- 
tions at random. 
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The standard test for superconducting coherence of a 
model superconductor has been to sample the helicity 
modulus 1/e: 

1 8tt 2 . 

Kk) = 1 -¥TT3 {qakqa - k) (5) 

In the limit k — > the phase transition is detected 
in the following way. For temperatures in the supercon- 
ducting phase 1/e 7^ and above the transition 1/e — 

In this letter we use an alternative test for the super- 
conducting transition namely the vanishing of the resis- 
tance |TT|| . The dissipation in a 3 dimensional super- 
conductor is caused by the creation of vortex loops and 
expanding them out to the system boundary. Allterna- 
tively if there is an external magnetic field that gives 
vortex lines through the system, the movement of these 
vortex lines will dissipate energy. The linear resistivity 
is defined by p = E/j for j — ► 0, where j is the ap- 
plied supercurrent density and E is the resulting induced 
electric field. The resistance R is given by the Nyquist 
formula 

R=^ dt{V{t)V{Q)). (6) 

The integral is evaluated as a sum over discrete time 
steps, defined as one MC trial move. The voltage V x (t) 
is defined by the fluctuation of loops and is calculated by 
the following procedure. N X+ (N X ^) denotes the number 
of accepted trial moves with a vortex-loop oriented in 
the x-direction as x + (x— ) for a MC sweep through the 
lattice. The + and — keep track of whether the vortex- 
loop is positively or negatively oriented. The voltage 
V x (t) at time t, in the ^-direction, is V x (t) oc N x+ —N x _. 
There are three resistances R x , R y and R z which all are 
equal in the isotropic case considered here. 

We consider now the finite size scaling. In three di- 
mensions 1/e obeys the scaling relation p3| , p4| ] 

L — ; s— r- as constant at T = T c and d — 3, (7) 

A finite size scaling relation for the resistivity, can be 
derived in the following way. The Josephson relation 

V ~ ±V* (8) 

relates the voltage to the time derivative of the gradient 
of the phase 6 of the superconducting order parameter 
p5f . From eq || we conclude that as T c is approached the 
voltage scales as V ~ 1/r where r is the dynamical time 
scale. Dimensional analysis of eq. |] leads to R ~ 1/r, 
where r is related to the correlation length through r ~ 
£ 2 . At T c the correlation length is cut off by the finite 
size L of the system and we have 

R ~ t- 1 ~ r 2 = L- z (9) 



In three dimensions we have the following relation for 
the resistivity, p = RL. Hence, the following finite size 
scaling relation for the resistivity: 

pL z ~ x « constant at T — T c and d = 3, (10) 

The Meptropolis algorithm does not in itself contain any 
reference to time. One can however show j2f| that there 
is a linear relation between time the scale of Langevin 
dynamics and Metropolis MC trial moves. The success 
of this similarity has proven itself in many simulations 

Now we turn to the results. The analysis is based on 
the finite size scaling relation eq. [l^. The temperature is 
measured in units of J x . The determination of z is done 
by the following minimisation procedure on our Monte- 
Carlo data. For a given value of z, we form the data 
curves p(L,T)i 2_1 as a function of temperature. Wc 
calculate the average separation St and S p between the 
crossings of these curves. The index T and p indicate the 
respective coordinates. For n curves there will be * 
crossings. The minimum S indicates the z for which the 
scaling relation eq. |l0| is full filled, and it determines 
the critical temperature T c . In figure [l] the functions St 
and S p are plotted versus z—1. The lattice sizes in the 
figure are L = 8, 10, 12, 14, 16. Both functions have a 
clear minimum, which occurs at nearly the same value 
z — 1 = 0.51. Less well converged data will not have co- 
inciding minima for the St and S p functions. We have 
also tried to exclude some of the lattice sizes in the cal- 
culation of St and S p but this does not change the result 
for z, at maximum 3% . Including lattice sizes L = 4 and 
6 will change the determination of z. Especially L = 4 
is outside the scaling regime and including both L = 4 
and 6 would change z to 1.51. The critical temperature 
is determined as the average intersection at the z that 
minimised St and S p and is found to be T c = 5.99. 

One might also note that if the data had not been well 
converged. The minimum in figure would have been 
less well pronounced. This is because the scaling expo- 
nent z — 1 is found to be small. For high temperatures 
there will always be the trivial scaling as there are no 
finite size effects in p for temperatures far above T c , and 
eventually one would find z = 1 far above T c . The in- 
set of figure [l] shows the calculated critical temperature 
versus the scaling exponent z — 1. From the inset we see 
that a large change of the scaling exponent z — 1 only 
gives a moderate change in T c . Taken together with the 
well defined minimum in St and S p we infer that the 
procedure to determine z is stable. 

In figure [l]b the finite size scaling is shown for 1/e in 
accordance with eq. ^. The evaluated critical tempera- 
ture corroborates the result achieved with the resistivity 
scaling. The critical temperature determined is in good 
agreement with determinations for the 3 dimensional Vil- 
lain model S . There are no adjustable parameters in this 



2 



procedure and we can clearly see there is a small finite 
size effect, as the the curves for larger lattices intersect 
at slightly lower temperatures. One might also note that 
as the scaling relation for 1 /e works it indicates that the 
statitic scaling exponents are the same as for the 3 di- 
mensional XY-mode\. 

In figure || the resistance scaling is shown for the z 
that minimised the spread in figu re [D . The data shows 
a very good splay at T c and eq. |10| is obeyed to high 
precision. The inset shows the resistivity as a function of 
temperature. From figure ^| it is evident that there is a 
small finite size effect. The curves for larger lattices cross 
at higher temperatures. The effect is small and T c will 
have its upper bound given from the 1/e scaling shown 
in figure [l]b. From the inset in figure |l|a an approximate 
value for z would be 1.5. 

We have used the Nyquist relation to determine T c di- 
rectly from the vanishing resistivity. From the size scaling 
near T c we determine the dynamical critical exponent z 
to be z = 1.51 ± 0.03. This result is interesting since it 
is equivalent to superdiffusive behaviour. Most models 
have subdiffusive behaviour, i.e. z > 2 pTj. It is also 
worth to emphasis that the result establish that the 3d 
XT-model and the 3d London lattice model has the same 
dynamical critical behaviour not only the same equilib- 
rium exponents. It is interesting to compare our result 
to a recent work by Wengel and Young |L6| a study of 
the Lattice superconductor in the limit A — was pre- 
sented. In this limit of the model they find z = 3. The 
difference between our result and their result is consis- 
tent with the fact that LLM for A < oo is not equivalent 
to the XT-model. 
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FIG. 1. Monte Carlo results for the LLM. Shown in figure 
a are results for the scaling relation eq. [l^. The functions 
S p (dashed curve) and St (solid curve) are drawn as function 
of the dynamical exponent z. The lattice sizes employed in 
the determination are L = 8, 10, 12, 14 and 16. The minimum 
occurs at z — 1.51. Lhe critical temperature of the system 
is determined to T c = 5.99. Lhe inset shows the determined 
T c as a function of z. In figure b the results for the scaling 
relation eq. [?] are shown. Lattice sizes L are 4 = stars, 6 = 
open circles, 8 = filled circles, 10 = open squares, 12 = filled 
squares, 14 = triangles and 16 = plusses. One can clearly see 
the curves for larger lattices intersect at lower temperatures. 



FIG. 2. Monte Carlo results for the scaled resistivity. The 
function p(T)L z_1 is plotted against temperature. The dy- 
namical exponent is determined from figure |l|a z = 1.51. Lat- 
tice sizes L are 4 = stars, 6 = open circles, 8 = filled circles, 10 
= open squares, 12 = filled squares, 14 = triangles and 16 = 
plusses. There is a finite size effect present, intersections for 
the larger lattices take place at a slightly higher temperature. 
The inset shows p as a function of T. 
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